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Abstract 
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■ We analyze the overlap of local color-octet meson operators with the T and the r\b 
CN . and their low-lying excited states, especially the first radial excitations. Our analysis 

is based on NRQCD and includes all terms up to order v A . We use a variety of source 
and sink operators as a basis for the variational method, which enables us to clearly 
separate the mass eigenstates and hence to extract the desired amplitudes. The 
results show the usefulness of the variational method for determining couplings to 
excited hadronic states. 
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^ . The inner structure of hadrons is one of the most shrouded and difficult is- 



sues of modern physics. Due to the confining nature of strong interactions, 
perturbative QCD calculations cannot describe physics at hadronic energy 
scales from first principles. However, lattice QCD, as a nonperturbative ap- 
proach, has succeeded in confirming and predicting properties of matter at 
these scales. Since a priori we are not well versed in the structure of hadronic 
states, we try to probe them on the lattice with suitable operators. In quan- 
tum field theory, the only condition for non-zero overlap of an eigenstate of 
the underlying Hamiltonian with a specific operator is that the operator has 
the same quantum numbers as the mass eigenstate. So all physical states, the 
ground state as well as all its excitations, are expected to couple non-trivialiy 
to currents which project out the corresponding quantum numbers. A recent 
paper by Liu and Luo [1], addressed to spectroscopy of charmonia, suggests 
that there is an exception to this rule. In correlators, which they constructed 
from local color-octet (or hybrid) operators, no low-lying excitations are vis- 
ible. Only the ground state and a much higher lying state, most probably a 
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hybrid excitation, appear. We want to check this result for bottomonium sys- 
tems in the pseudoscalar and vector channel. Furthermore, we want to provide 
quantitative estimates for the couplings (another approach was recently used 
for excited pions [2]). To do so, we use the variational method, which is based 
on the construction of a cross correlator matrix. This approach enables us not 
only to investigate the spectrum, but more importantly, to obtain ratios of 
couplings of different local operators to a physical state. 

Of course, none of these operators, which merely serve as an expanding basis 
for the states, are physical in QCD, since mixing between them is possible. 
In the framework of NRQCD we can actually identify the terms in the La- 
grangian, which are responsible for this, so called, configuration mixing. This 
has been studied in [4], for instance. Since such configuration mixings are 
ever-present in physical hadronic systems, all states should have overlap with 
operators which have the same quantum numbers, regardless of whether they 
are hybrid in nature. Since hybrid operators contain elements of the field 
strength tensor, the couplings associated with them are divergent quantities. 
Rather than removing these divergences, it is our aim here to show that, at 
typical lattice scales, the hybrid couplings are finite, yet relatively small. 

This paper is organized as follows. Section 2 describes the variational method 
we use for extracting masses and couplings. Our implementation of the NRQCD 
framework to propagate the fermions is discussed in Section 3. Section 4 gives 
an overview of the configurations on which we are running our simulation. 
The actual results are presented in Section 5 and Section 6 contains conclud- 
ing remarks. 

A preliminary report of these results can be found in [3]. 



2 Method 

To gain information on the couplings we rely on the variational method. The 
starting point is the construction of a cross correlator matrix CV,-: 



where X n is the eigenvalue corresponding to the eigenvector if; 71 . Since, due to 
fluctuations, is not exactly symmetric (although it is within errors), we 



C y (t) = (0|O i (t)O i (0)|0>. 




Then we can solve the generalized eigenvalue problem 




(2) 
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symmetrize it by hand in order to make the diagonalization procedure more 
stable. 



The eigenvalues are then given by [6,7] 



\ n (t,t ) = A n e- E ^-^[l + O(e- AE ^-^)}, (3) 

where E n denotes the mass of the nth state and AE n the mass difference to 
the state closest in mass to the state n. This correction is due to the use of a 
finite number of operators in the basis. However, for large enough values of t, 
we have a single mass eigenstate in each channel. Thus, the variational method 
enables us to clearly separate the ground state, many lower lying excitations, 
and even ghosts [8] (which do not play any role in this work due to the large 
mass of the quarks). 

For our purposes here, even more important than the eigenvalues is the analy- 
sis of the eigenvectors. As a starting point consider the spectral decomposition 
of C: 

CM = E ^~ Ent (4) 

n=l ltj n 

where i>f = (0|Oj|n) = (Oi\n), that is the overlap of ith operator with the nth 
state. To extract this matrix element consider the right-multiplication of the 
correlator matrix with one of the eigenvectors: 

E Ca (t) r 3 = ^r<a n e- E "*[l + 0(e" A ^)] *** ^aV w , (5) 

where a n = Y,iV^*ip" '. For large enough times a single mass eigenstate is 
projected out of the sum (see Eq. (3)). We can now left-multiply with the 
same eigenvector to obtain 

E#*^)^^K|V^. (6) 



n 



In the event that the correlator matrix is purely real (or purely imaginary) we 
can calculate the a n and divide them out from Eq. (5) to obtain the couplings 
i>". Since we are only interested in ratios of couplings, we proceed with Eq.(5). 
To cancel the operator independent amplitude, the energy prefactor and the 
exponential we construct a ratio of the following kind: 

XjCijM? t>>t0 (O t \n) = Zr^Oiln)^ 

EiC kl (t)W (O k \n) Z k l (O k \n) r J 1 ) 

where Zj is the renormalization constant for the zth operator and (. . .) 
denotes the renormalized quantity. 
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It is important to note that in using a finite, non-orthogonal basis of operators, 
the individual overlaps may not be well defined. It is vital that one checks, for 
example, that the above ratio of local couplings does not change when further 
operators are added to the basis. 

The success of the variational method strongly depends on the choice of op- 
erators one includes in the basis. Unfortunately, there will in any case only 
be a finite number of operators available, which can just span a subspace of 
the Hilbert space of states. Therefore, one tries to use operators which have 
small relative overlap in order to explore a large area of the space the me- 
son wavefunctions live in. Since we want to investigate the rjb and the T, we 
use pseudoscalar and vector currents. For both we have a "normal" and a 
"hybrid" version. Table 1 gives an overview of the local operators we use. It 
additionally contains operators coupling to states which we only address in 
the spectrum section. The P-wave states are moreover needed to set the scale. 
In order to assemble our basis with more linearly independent operators, we 
additionally smear the quark and the antiquark field independently with two 
different smearing levels n naTrow = 16 and n wi d e = 48 using gauge invariant 
Jacobi smearing with k = 0.2. 

In total, twelve different operators, each at the source and the sink, are avail- 
able for constructing the cross correlator matrix. 



Table 



state 




normal operator 


hybrid operator 


Vb 


o-+ 


xV 


X^iBi(f> 


T 


1~ 


x'W 


X ] B l <f ) 


Xbo 


0++ 


x^i D i4> 




Xbi 


1++ 


X^ijk<?jD k (j) 




Xb2 


2++ 


xH^iDj + ajDi - l6ij<T k D k )<j> 




h b 


1+- 


X f A0 




D-wave 


2-+ 


xH D i D j + DjDi - l5ijD k D k )(p 




exotic 


1-+ 




X^ijk<?jB k (f> 


1 



Overview of the used operators Oj . The meaning of <\> and x will be defined in the 
subsequent section. 
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3 NRQCD 



The calculation of the propagators is performed in the framework of NRQCD 
(see [9]), which is perfectly suitable for our bottomonium systems, where the 
quarks move with small relative velocities. We include all terms up to 0(v A ) in 
our NRQCD Lagrangian, where v is the velocity of a quark, according to the 
power counting in [9]. Since we are working in a nonrelativistic approxima- 
tion, it is convenient to separate the Dirac four spinor ip into two sets of Pauli 
spinors, for the quark and \ f° r the antiquark. Furthermore, the propaga- 
tion of the fermions can be described by the quantum mechanical evolution 
operator for imaginary time e~ Ht . By expanding this operator we obtain for 
the propagation: 



0(y, t + «) = (l - a -M±^j " Ut(t) (l - 2^) " (1 - aSHtt)) «x, t), 

(8) 

where H is 



H 



A a A 2 
2M ~ 4^4M 2 



and SH is 



(9) 



5H = 



8M 3 

w 

8M 2 
9 



(10) 



(V-E-E-V)-^.(VxE-ExV) 



2M 



a B. 



The tildes denote improved versions of the corresponding derivatives. We use 
n — 2, which is more than sufficient in our case. B and E are the magnetic and 
electric fields created via the usual clover formulation. The last two terms of 
(10) are responsible for the configuration mixing mentioned above; however, 
it turns out that the a ■ B term produces a much larger effect. 



To correct for tadpole contributions we divide each link by a factor u which 
is given by the fourth root of the average plaquette [10] . The quark mass for 
our simulation is determined from finite momentum correlators for the T by 
tuning the kinetic mass extracted from the non-relativistic energy-momentum 
dependence to the experimental mass of the T. 
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4 Configurations 

We are working on configurations provided by the MILC-collaboration [11]. 
They were generated using improved staggered fermions and the Liischer- 
Weisz gauge action. Table 2 shows the parameters of the lattices used. For 
the lattice spacing, there are two values given. The first one comes from the 
analysis of the spin-averaged T IP-IS splitting, the second one is given by 
the MILC-collaboration, where they used the improved Sommer parameter r\ 
to set the scale. Note that for the quenched lattice we are slightly above the 
values from the MILC-Collaboration. This is most likely due to the wrong 
curvature of the qq potential in quenched simulations, thus influencing the P- 
and the S-waves differently. 

To obtain the physical fe-quark mass at the corresponding scales, we inter- 
/extrapolate to the experimental T mass. However, both for the spectra and 
the ratios of the couplings, we find a very weak dependence upon the quark 
mass parameter. 





o _1 [MeV] 


volume 


N f 


Q>rn sea 


arrib 


m p b hys [MeV] 


# of configs. 


8.40 


2378(7)/2279 


28 3 x 96 





OO 


1.7, 1.8 


4114(72) 


160 


7.09 


2097(9)/2252 


28 3 x 96 


2+1 


0.0062/0.031 


1.7, 1.8 


4026(65) 


210 


6.76 


1495(6)/1587 


20 3 x 64 


2+1 


0.01/0.05 


2.4, 2.5 


4111(48) 


410 



Table 2 



The three different lattices we use in our simulations together with the correspond- 
ing parameters. The first value for the inverse lattice spacing comes from the IP-IS 
splitting, the second one is given in [11]. 



5 Results 

5. 1 Spectra 

At first we want to present the results for the masses in the considered chan- 
nels. Several efforts have already been made to investigate the spectrum of 
bottomonia ([12], [13], [14], [15]). Table 3 gives some of the S-wave splittings 
for the fine quenched and dynamical lattices. Figure 1 shows the plot of the 
absolute masses. The fits are one exponential fits to the corresponding cor- 
relators of the states, using the full covariance matrix. For the pseudoscalar 
and the vector channel we use the variational method in the basis Nil, Nln, 
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Nnn, Nww, Hll and thus are able to obtain masses of their excited states 
reliably, the masses of the other states are extracted from single correlators. 
The absolute mass offset is fixed by the value of the T. For both types of 
configurations the quark mass parameter was set to amj = 1.7. It is worth 
noting that the experimental input for the r)b comes from a single event, so it 
is not very reliable. 
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Fig. 1. Spectra of the available states, obtained by fitting the corresponding corre- 
lators, for both quenched(/3 = 8.40) and dynamical configurations (/? = 7.09). 



The capital letter denotes the type of the operator: N=normal, H=hybrid; the two 
small ones give the smearing of the quark and the antiquark respectively: l=local, 
n=narrow, w=wide. 
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Fig. 2. Effective mass plots of the available states from the quenched lattice with 
(3 = 8.40. The results for ^ and 1 come from the variational method, the other 
effective masses have been obtained from single correlators. The blue lines indicate 
the fitted plateaus. The fitting errors can be told from their width. Note the different 
scale for the exotic 1 . 
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Fig. 3. Same as in Fig. 2 but for the dynamical lattice with (3 = 7.09. 
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0-+(15) 


0-+(2S) 


0" + (3S) 


1--(1S) 


1-(2S) 


1~(3S) 


experimental 


9.300(28) 


- 


- 


9.460 


10.023 


10.355 


quenched 


9.389(4) 


10.004(10) 


10.563(39) 


9.460(4) 


10.031(7) 


10.487(16) 


X 2 /dof 


2.49/5 


7.42/4 


2.50/4 


6.76/5 


2.36/4 


5.64/4 


fit range 


4-10 


3-8 


3-8 


4-10 


3-8 


2-7 


dynamical 


9.423(5) 


9.970(11) 


10.485(82) 


9.460(6) 


9.983(15) 


10.409(78) 


X 2 /dof 


4.56/5 


6.56/5 


2.43/4 


2.65/4 


1.71/3 


1.27/3 


fit range 


4-10 


4-10 


4-9 


10-15 


7-11 


5-9 




0++ 


1++ 


2++ 


1+- 


2-+ 


1-+ 


experimental 


9.859 


9.893 


9.912 








quenched 


9.878(5) 


9.903(5) 


9.927(7) 


9.956(5) 


10.308(7) 


11.385(268) 


X 2 /dof 


2.79/5 


4.52/4 


2.44/2 


5.89/5 


8.03/7 


0.74/2 


fit range 


10-17 


12-18 


16-19 


10-16 


7-15 


3-7 


dynamical 


9.860(8) 


9.884(8) 


9.907(7) 


9.891(8) 


10.214(17) 




X 2 /dof 


3.15/4 


5.09/4 


3.04/4 


3.02/3 


4.63/4 




fit range 


11-16 


11-16 


11-6 


11-15 


9-14 





Table 3 



Fitting results for the masses in GeV for both quenched(/3 = 8.40) and dynamical 
configurations (/3 = 7.09). 
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5.2 Starting basis 



Local operators are well defined in the context of quantum field theory and 
therefore the ratio of their couplings is of particular interest. Furthermore, this 
will provide us with a clear answer if radial" 2 ! excitations have overlap with 
local hybrid operators. To access these ratios we proceed to the analysis of the 
eigenvectors of the cross correlator matrix. As described above we construct 
the correlator matrix with a variety of source and sink operators and diagonal- 
ize it. By this procedure we hope to clearly separate the mass eigenstates and 
use the eigenvectors to shed light on their overlap with the included operators. 

We start our analysis for the T, which is probably the experimentally more 
interesting state, on the dynamical lattice with (3 = 7.09 and a quark mass 
parameter of ami, = 1.7. Later on we will also investigate the other lattices 
with different parameters to see what changes arise. 

The smallest basis, including a hybrid operator, which reveals reasonable re- 
sults, is Nll(l), Nnn(2), Hll(3). The numbers in brackets are given just for the 
sake of clarity in the coupling ratio plots. The effective masses of the eigen- 
states are shown in Fig. 4. The ground state and the first radial excitation are 
clearly visible. 

To check if our diagonalization worked correctly, we can reconstruct the eigen- 
value by multiplying the eigenvector with the cross correlator matrix Cij(f)i]f-. 
The indices i,j label the operators and n the states. This product should 
have the same exponential behavior, i. e., the same mass, as the corresponding 
eigenvalue. Of course, they may differ in their amplitude since Cij(t)vp™ gives 
the eigenvalue times some overlap factor. But much more important is the 
fact that these overlap factors ultimately provide us with the ratios of the 
amplitudes. To give a quantitative estimate for the ratio, we have to plot the 
ratio with respect to time. When the ratio has plateaued, we can be sure that 
higher excited states do not play a role any more. 

Figure 5 shows the ratio of the overlaps of the local hybrid and the local nor- 
mal operator for the ground state and the first radial excitation. We clearly 
find that radial excitations have non-zero overlap with local hybrid operators. 
This is obvious since the radial excitations have the same quantum numbers 
the hybrid operator couples to. In both cases the coupling of the local nor- 
mal operator to the state is about 90 times larger than the coupling of the 
local hybrid operator. Of course, these numbers still require renormalization. 

2 The low-lying excited states of bottomonia of definite J PC are "mostly" radial 
excitations. Of course, configuration mixings will ensure that some small admixtures 
of other configurations (gluonic, orbital, etc.) are present. However, in the following 
we use "radially excited" and "excited" interchangeably. 
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However, this is beyond the scope of this work since we are only interested in 
showing that the couplings to hybrid operators are finite, yet small for typical 
lattice scales, and that the variational method can be used to extract these 
couplings. By constructing a ratio of these ratios we define the R value: 

pnn' = (°i\ n )ren { i\ n ')ren /i i \ 

lk (O k \n) re J (O k \n') re J [ } 

where, of course, % and k are not summed over. According to Eq. (7) R can 
be expressed in terms of the cross correlator matrix and its eigenvectors: 

The renormalization constants which only depend on the current cancel and 
a further cutoff dependent factor which may come from disagreeing mass di- 
mensions of the operators cancels too. If one of the operators is a hybrid one, 
one would also expect a scale dependence due to its nonultralocal nature, 
more precisely F^ v appearing in the hybrid operator is constructed from a 
2x2 clover. We assume the factorization of this dependence and thus Eq. (11) 
provides a scale independent physical quantity. We will actually show this in 
Section 5.4. Therefore, in order to remove any cutoff dependence, we utilize 
Eq. (11) and obtain as a physical result that the ratio is about equal for the 
ground and the first excited state (see Fig. 6). The fit to a constant from 
t = 3 to t = 10 yields B%{ = 1.0464(63), where the error comes from single 
elimination jackknifing. 
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Fig. 4. Effective masses of the three eigenvalues in the basis Nll(l), Nnn(2), Hll(3) 
for the dynamical lattice with (3 = 7.09 and ami, = 1.7. to is the normalization 
timeslice. 
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P=7.09, 28 x96, N f =2+1, am b =1.7, t =l 
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Fig. 5. The ratio of the couplings of the two local operators to the ground state is 
shown in the left plot, the same for the first excitation is in the right plot. 
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5.3 Enlarged basis 



In order to obtain better signals from our correlators and to show a general 
feature of our operators, we increase the number of normal operators in the 
basis. The enlarged basis is built up from Nll(l), Nln(2), Nnn(3), Nww(4), 
Hll(5). Figure 7 shows the effective masses of the eigenvalues. For this choice 
of basis we manage to obtain an acceptable signal even for the fourth state. 
The ratios of the couplings for the ground state and the first radial excitation 
Fig. 8. Figure 9 shows the ratios for the third and fourth state, two further 
excitations which have been skipped in the previous basis. Again the results 
are close to those for the ground state and the first radial excitation. The 
value of R evaluated in this basis for the ground and first excited states is 
shown in Fig. 10. The fit to a constant in the range t — 3 to t — 10 gives 
Rlf = 1.0348(87), confirming the result from the previous basis. 

A property of local operators reveals itself when we look again at the ratio of 
amplitudes. A comparison of the ground state and the first radial excitation 
ratio in two different bases is shown in Fig. 11. The explicit forms of the 
bases are: Basis 1 - Nll(l), Nln(2), Nnn(3), Nww(4), Hll(5), basis 2 - Nll(l), 
Nln(2), Nnn(3), Hll(4), Hln(5), Hnn(6). Even though, we change our basis by 
substituting smeared normal operators by other smeared ones, the ratios of 
the couplings of the local operators remain the same for these states. This 
implies that local operators are "approximately" orthogonal to smeared ones 
(i. e., they have very small spatial overlap). If they would overlap significantly, 
adding a smeared operator to the basis would change the ratio, because the 
smeared operator could "steal" some contribution of the projection of the 
physical state onto the local operator. When looking at ratios of couplings 
of two smeared operators, we see that they change significantly when another 
smeared operator is added to, or removed from, the basis. That means that the 
smeared operators we use are not orthogonal to each other. Since the coupling 
ratio of the local operators stays the same in every arbitrary basis, as we found, 
we can also exclude that the smeared operators occasionally steal exactly the 
same contribution from the local normal and the local hybrid operator in such 
a way that the ratio remains the same. The result about the orthogonality 
of local and smeared operators is very plausible. The spatial width of the 
local operators is near zero and they have finite height. So their convolution 
with the Jacobi smeared operators should be quite small, as we found. The 
constancy of the ratio under the change of basis is an important justification 
for our conclusions. 

It is worth noting that including further smeared operators in the basis does 
not help to improve the outcome; quite the contrary, they enhance the overlap 
with higher excited states or they contribute more noise than new information 
and thereby disrupt the signals. 
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Fig. 7. Effective masses of the five eigenvalues in the basis Nll(l), Nln(2), Nnn(3), 
Nww(4), Hll(5) for the dynamical lattice with (3 = 7.09 and am& = 1.7. 



(5=7.09, 28 x96, N f =2+1, am b =1.7, t Q =l 



0,02 
0,018 - 
0,016 - 
0,014 - 
0,012 - 

0,01 - O 
0,008 - 
0,006 - 
0,004 - 
0,002 - 



I 1 I 1 I 



I ' I 







I i I i I 



I 1 I 1 I 



I i I i I 



2 4 



8 10 12 14 2 4 6 8 10 12 14 



0,02 
0,018 
0,016 
0,014 
0,012 
0,01 
-£f0,008 
0,006 
0,004 
0,002 




Fig. 8. The ratio of the couplings of the two local operators to the ground state and 
to the first radial excitation. Both in the new 5x5 basis. 
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Fig. 9. The ratio of the couplings of the two local operators to the third state is 
shown in the left plot, the one for the fourth state is in the right plot. Both in the 
new 5x5 basis. 
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P=7.09, 28 x96, N =2+1, am =1.7, t =1 
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Fig. 11. The ratio of the couplings of the two local operators to the ground state 
is shown in the left plot, the one for the first excitation is in the right plot, both 
for two different sets of operators [basis 1 is Nll(l), Nln(2), Nnn(3), Nww(4), Hll(5) 
and basis 2 Nll(l), Nln(2), Nnn(3), Hll(4), Hln(5), Hnn(6)]. The datapoints for the 
second basis are shifted by 0.3 for the sake of clarity. 
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5.4 Scale dependence/ quenching effects 



To verify that R defined in Eq. (11) really is scale independent, as we claim, 
we do the same calculation in the previous 5x5 basis for a coarser lattice 
(a -1 = 1497 MeV compared to a -1 = 2101 MeV in the previous calculation) 
with arrib = 2.4. The results are shown in Fig. 12 and Fig. 13. From Fig. 12 
follows a clear scale dependence for the individual couplings. However, Fig. 13 

1 2 

certifies that this dependence vanishes for R. There R 5 \ fitted to a constant 
in time at three different scales, two from dynamical and one from quenched 
lattices, is plotted. Just considering the case with sea quarks included we see 
that the two points are consistent within the errors, proving that any scale 
dependence is very small. To guide the eye we performed a fit to a constant 
with the two dynamical points, yielding 1.00784(836), which is consistent with 
one. The fact that the quenched result is relatively low shows the different 
effect of virtual quark-antiquark pairs for the ground and the first excited 
state. However, the hybrid coupling is also seen to be small on the quenched 
lattices. 
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Fig. 12. The ratio of the couplings of the two local operators to the ground state is 
shown in the left plot, the one for the first excitation in the right plot. Each plot is 
for two different lattices with a -1 = 2101 MeV (fine) and a -1 = 1497 MeV (coarse). 
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Fig. 13. Plot of R%1. The green line shows the result for fitting the dynamical points 
to a constant. The fitting error can be told from its width. 
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5.5 Further results 



This section briefly covers the results for the 775, which is, from an experimental 
viewpoint, less accessible, since its ground state has been observed in only one 
type of experiment and none of its excitations have been detected up to now. 

Figure 14 shows the effective masses of the eigenvalues in the basis Nll(l), 
Nln(2), Nnn(3), Nww(4), Hll(5). The local coupling ratios are plotted in Fig. 
15. Typically, the signals for pseudoscalars are slightly better than the ones 
for vectors. The increase of the ratio of the local couplings by about a factor 
three can be traced back to the fact that in the h hybrid operator all three 
components of the .B-fleld are included. 
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Fig. 14. Effective masses of the five eigenvalues in the basis Nll(l), Nln(2), Nnn(3), 
Nww(4), Hll(5) for the pseudoscalar on the dynamical lattice with (3 = 7.09 and 
arrib = 1.7. 
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Fig. 15. The ratio of the couplings of the two local operators to the ground state is 
shown in the left plot, the one to the first excitation in the right plot. Both for the 
pseudoscalar state on the dynamical lattices. Additionally the results for fitting the 
ratios to a constant with and without diagonalization are shown. The fit range for 
the direct fit of the (3 x 3) correlator matrix is t = 13 — 18 for both states. 
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As a crosscheck we also determined the couplings of the two local operators to 
the two lowest lying states by directly fitting a three by three correlator matrix. 
To do so, we build up a matrix of correlators with three different operators, 
namely Nll(l), Nnn(2), Hll(3). Then we fit the matrix to the functional form 

C t ,(t)=A^A J e- aE ^ + AlA]e- aE ^ t , ij = 1,2,3, (13) 

where A® is the coupling of the ground state to the ith operator and Aj is the 
coupling of the first excited state to the zth operator. We find that fitting the 
above form with a larger basis of operators and with more masses becomes 
highly problematic. 

The results for the ratios of couplings and the corresponding R value obtained 
by direct fitting together with the ones after diagonalization in the basis Nll(l), 
Nln(2), Nnn(3), Nww(4), Hll(5) are given in Table 4, both for the vector and 
the pseudoscalar channel (see also Fig. 15). The time fit range to determine 
the couplings in the case of direct fitting is t — 13 to 18 for all three quantities, 
in the case of the fit after diagonalization t — 6 to 13 for the ground state and 
t — 5 to 13 for the excited state and R. The range for direct fitting multiplied 
by the nine different operator combinations yields 54 degrees of freedom to 
determine eight parameters, namely six couplings and two masses. The errors 
in both cases were calculated by single elimination jackknifing. The numbers 
suggest that the direct fitting method is slightly better for the ground state, 
but much worse for the excited one. One may criticize such a comparison 
by pointing out that the direct fitting takes place with only three operators 
in the basis, while the diagonalization works with five. But actually this is 
the main virtue of the variational approach. With the diagonalization, one is 
able to work with more operators, thus putting in more information. In turn, 
this provides eigenvalues which more readily correlate single mass states. One 
can then start fitting (single exponentials) much earlier, thereby avoiding the 
noisier areas of the time range. However, we note that the outcomes of both 
methods show reasonable agreement. 



6 Conclusions 



In the end of this work we want to summarize our results and draw conclusions 
from them. We found, as do many others working on spectroscopy, that the 
variational method is a very promising approach since it allows one to clearly 
separate the individual mass eigenstates of the theory and, furthermore, to 
find out the ratio of the couplings of a state to different operators included 
in the basis of the cross correlator matrix. The low cost of simulations in the 
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pseudoscalar 


vector 


direct fitting 


diagonalization 


direct fitting 


diagonalization 


A V A i 

A l(5)/ A l 

R 


0.033628(50) 
0.03524(231) 
0.9634(740) 


0.033416(122) 
0.03313(47) 
1.0082(99) 


0.011381(23) 
0.011993(728) 
0.9519(608) 


0.011372(48) 
0.011024(142) 
1.0348(87) 



Table 4 



Results for the direct fit of the correlator matrix in the basis Nll(l), Nnn(2), Hll(3) 
and for the fit after diagonalizing the correlator matrix in the basis Nll(l), Nln(2), 
Nnn(3), Nww(4), Hll(5). For the direct fit we state R.l\, whereas for the fit after 

1 2 

diagonalization R 5 \. 



NRQCD framework also contributes to the success of our analysis. 

We find the quite amazing fact that the ground state and the first few radial 
excitations of bottomonia have about the same ratio for their couplings to 
the two local operators we considered. The outcome for the radial excitations 
clearly contradicts the results of [1]. They cannot find the first radial excita- 
tion in their correlator built up from only hybrid currents, not because it is 
"skipped," but because the coupling between this state and the hybrid oper- 
ator is too small. With inputs of single symmetric correlator runs we actually 
are able to conclude that in the correlator built from hybrid operators the 
amplitude of the radial excitation is down by 0(100) compared to the hybrid 
excitation. The smallness of the coupling is the only reason why one is not 
able to see any radial excitation in hybrid correlators. 

A very important condition for the functionality of our approach is that the 
ratio of the couplings of the two local operators is independent from the choice 
of basis. In other words, the local operators have to be orthogonal to the 
smeared ones used in the basis. We also show that the R value, on which we 
mainly base our statements, is rather independent of the scale. However, slight 
sea quark effects can be observed for R. 

As a crosscheck for our results we directly fit the correlator matrix without any 
diagonalization. Although the quality of the fit is comparable for the ground 
state, it is inferior for the first radial excitation. Separating the individual 
eigenstates by diagonalization enables us to start the fit much earlier, about 
t of 5, whereas for the direct fitting the starting time is about t = 12, leading 
to larger errors for the radial excitation. 
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